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We study the behaviour of Wilson and 't Hooft loop operators for the 2+1 dim. Abelian-Higgs 
model with Chern-Simons term. The phase of topological symmetry breaking where the vortex field 
condenses, found by Samuel for the model in the absence of the Chern-Simons term, persists in 
its presence. In this phase, the topological linking of instantons, which are configurations of closed 
vortex loops, with the Wilson loop on one hand and with the 't Hooft loop on the other hand 
ON , gives rise to a long-range, logarithmic, confining potential between electric charges and magnetic 

flux tubes. This is surprising since the gauge field is short range due to the Chern-Simons term. 
Gauss' law forces the concomitance of charge and magnetic flux, hence the confinement is actually 
of anyons. 
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I. INTRODUCTION 



The Abelian-Higgs model has long been a model of wide interest from condensed matter to particle physics. 
It is the effective theory for the description of superconductivity with a charged condensate. It also provides the 
simplest example of the Higgs' mechanism where spontaneous breaking of a gauged symmetry evades the Goldstone 
theorem and provides for a purely massive spectrum. The theory is expected to have two phases, the Higgs (or 
superconducting) symmetry broken phase and the Coulomb phase. The expectation value of the Wilson loop operator 
which measures the potential between external electrically charged particles, distinguishes between the two phases. 
In lower dimensions it is possible that non-perturbative phenomena disturb this straightforward conclusion. In 1+1 
dimensions instanton effects actually conspire to restore the symmetry in the putative Higgs phase, and indeed there 
reappears only a Coulomb phase with however a distinctly different strength for the linear confining potential [|l). 
The result is essentially topological, the Wilson loop divides space-time into two disjoint parts, the Wilson loo_ 



loop 



integral measures whether the instanton is (topologically) inside the loop or outside the loop. In 2+1 dimensions 
• . nonperturbative instanton contributions to the expectation value of the Wilson loop were computed. In the Higgs 
region topological solitons exist in the spatial plane as vortices, their continuation to three dimensional Euclidean 
space are vortex loops which serve as the instantons. Their contribution to the path integral requires summing over 
all configurations of vortex loops. Again the contribution is essentially topological, the Wilson loop integral measures 
the topological linking number of the vortex loop with the Wilson loop. 

On a lattice of fixed spacing one finds the contribution of loops of length L to be approximately e~ mL u L , where 
m is the vortex mass per unit length, and n is an intrinsic parameter such that [i L counts the degeneracy of loops 
of length L. A phase transition occurs when m — In /i ~ where vortex loops of all lengths become important, 
which has been seen to occur in lattice simulations of the model, see e.g. [pi. A partition function of a gas of loops 
can be converted to a effective field theory of a scalar field called the vortex field. Condensation of the vortex field 
corresponds to the "spaghetti vacuum" and is called topological symmetry breakdown. The vacuum is populated 
of vortex/anti-vortex pairs restoring the U(l) symmetry giving back the Coulomb phase. In the region where the 
vortex field condenses the vortex loops induce a nonperturbative logarithmic confining force between external charged 
particles g). If the external particles have charge q, the strength of the logarithmic potential is periodic in q with 
period e, the charge of the elementary scalar particles in the theory, since they can screen the external particles when 
q is a multiple of e. 

In 2+1 dimensions or more, one can also look at the behaviour of the 't Hooft loop operator Q. It is defined by 
the action of a singular gauge transform on the fields, and at each time slice it creates a vortex/anti-vortex pair. The 
behaviour of the 't Hooft loop expectation value describes the inter-vortex potential. The perturbative and vortex 
loop configurations give a perimeter law dependence to the 't Hooft loop for both phases of the theory. 

If the Chern-Simons term (CS), with coefficient k is included we ask what is its effect on the Wilson and 't Hooft 



loop expectation values? The immediate problem that arises is that the CS term becomes imaginary in Euclidean 
space. Since the action is complex, the Euler-Lagrange equations over determine the system and there is in general 
no solution to the field equations (in terms of real fields). There do exist complex solutions however, for small k we 
expect that the real vortex loop configurations found at k — will continue to be the most relevant configurations. 
Our approach is thus the same as [E|, who considered the effect of the CS term to the monopole contribution in the 
Georgi-Glashow model. In fact we can show complex vortex solutions exist in the present case but we do not use 
them here. 

We find that for the 2+1 dimensional Abelian-Higgs model the addition of the CS term does not significantly alter 
the vortex loop contribution to the Wilson loop expectation value for small k. The behaviour of the 't Hooft loop 
is however significantly altered. In a manner similar to the case without the CS term, there is a phase transition 
at the point where the entropy factor of the vortex loops overcomes the Boltzmann factor. When the vortex field 
condenses, vortex loops give a logarithmic potential between both the external electric charges and vortices. This 
result is very surprising given that the photon is now massive and perturbatively no long range forces exist. It is 
unclear to us whether this behaviour persists throughout the unbroken symmetry phase. For the vortex entropy factor 
to overcome the Boltzmann factor the classical mass of the vortex must be smaller than some critical value where the 
phase transition occurs. But topological vortex configurations exist classically only for a certain range of parameters. 
The mass of the vortex can be continuously reduced to zero, in this limit the topological stability of the vortex is 
lost. This corresponds to the region where the gauge symmetry is unbroken classically. Since we cannot construct 
topological vortex configurations here, our methods do not extend to this region. We speculate that the confining 
logarithmic potential continues throughout the symmetry unbroken phase. 

In Section 2 we describe the model and its vortex solutions. We then describe the contributions of the vortex 
loops to the path integral and discuss how to treat the path integral with the inclusion of the complex CS term in 
Euclidean space. In Section 3 we compute the instanton contribution to the Wilson and 't Hooft loop expectation 
values and conclude with comments. 



II. THE ABELIAN-HIGGS MODEL WITH CS TERM 

The Abelian-Higgs model with CS term in (2+1) dimensions has Lagrangian density given by 

C = {- J V" + \d^ (D^y - ^(|0| 2 - r, 2 ) 2 + ^'V^} - (2- 1 ) 

here <f> is a complex scalar field, A^ is a U(l) gauge potential, D^ = d^ — ieA^, F M „ = d^A v — duA^ (/j,, v = 0, 1, 2) 
with A, e, k and rj 2 coupling constants, k, A, rf and e 2 each have the dimension of mass. We are considering 
noncompact U{\) gauge theory, there is no quantization condition on k. When k = 0, if rf < the gauge symmetry 
is unbroken and the perturbative spectrum consists of a massless photon and a massive charged scalar particle of 
charge e. If if > the gauge symmetry is spontaneously broken and the particle content is given by a massive vector 
particle and a massive neutral scalar particle. In addition classical time independent solitons exist in the form of 
Nielson-Olesen vortices. The vortices carry a topologically conserved integer valued charge usually called the vortex 
number. The vortex mass m is approximately equal to rj 2 . For rf < the topological vortex solution no longer exists. 

If the CS term is now included, k^O, the theory is altered in that the photon is massive even when rj 2 < 0. 
The CS term automatically gives the photon a mass proportional to k, this is called topologically massive QED. In 
the unbroken phase when rf < the perturbative particle spectrum consists of a massive charged scalar field and a 
massive gauge particle. When the U(l) symmetry breaking occurs, for rf > 0, the particle content now consists of a 
massive neutral scalar particle and two massive gauge particles. In this case vortex solutions are also present. Gauss 
law implies that any vortex solution must have electric charge Q (and vice versa): 

d 2 x J° = Q = k$ , (2.2) 

where $ = J d 2 xFn is the flux and J^ is the conserved electro-magnetic current. Both vortices and charged particles 
are anyons. The first solutions were proposed by Paul and Khare, |J, a radially symmetric charge one vortex is of 
the form 



^ = f(ry 9 , A(r) = -^±e e , A Q (r,6) = A (r) . (2.3) 

r 

f(r) and A(r) satisfy boundary conditions required by finite energy and Ao(r) satisfies Ao(Q) — and A(oo) = 0. 
f(r), A(r) and Aq(t) can be solved for numerically. 

The Wilson loop operator is defined by 



W{C) = exp(iq <b A t dx l ) , (2.4) 

where C is a closed oriented loop in Euclidean three-space, x l (s), and q a real number. The Wilson loop measures the 
potential between external electric charges of charge q. The 't Hooft loop operator M, B(C), is the dual object to the 
Wilson loop operator. Its expectation value measures the potential between external flux tube solitons, and especially 
the vortices. Its action on a field eigenstate |.Aj(x), </>(x) > transforms this state by a singular gauge transformation 
into \A t c (x), (f) Qc (x) >. The gauge transform VL c {x) with 

ftc(x) = exp(iu)c(%)) , (2.5) 

it is defined by the following. If another closed curve C" links with C n times in a certain direction and C corresponds 
to x l {6), < 9 < 2tt, then uc(x(2tt)) — uic(x(0)) + 2nn. The singular gauge transform creates a flux tube along 
the curve C, as measured by J c , A i c (x)dx l . Thus B(C) creates an infinitely thin bare flux tube along the curve 
C. To calculate its expectation value it is necessary sum over all field configurations which include this infinitely 
thin flux tube, and each configuration of action A is weighted with e . Is it not possible to vary the flux of the 
vortex/ anti-vortex pair (except in integer multiples) that is created by the 't Hooft loop, unlike the Wilson loop where 
the charge of the external particle q can be arbitrary. The flux of the bare vortices is always topologically quantized 
in integer multiples of 2tt /e. 

Perturbatively the expectation values of both W(C) and B(C) have perimeter law behaviour for large curves C. 
We take C to be rectangular of lengths t, and r, with t ^> r then their expectation values have the following behaviour, 

< W(C) >« e- v{r)t , < B(C) >« e - v ' (r)t . (2.6) 

V(r) and V'(r) measure respectively the potential between external charges and fluxes. The zero order approximation 
for V(r) can be shown to reproduce the classical electrostatic potential in two dimensions, which for large r behaves 
as Inr for a massless gauge theory and e~ mAr /(niAr) 1 / 2 for a massive gauge theory where tua is the gauge mass. 
This remains a good approximation for weak coupling, defined by e 2 <c M, where M is the smallest mass in the 
theory. For the 't Hooft loop, V'(r) rs E where E is the energy of the created object, in this case E ss 2m where m 
is the mass of the vortex. This implies the 't Hooft loop has perimeter law behaviour perturbatively recovering the 
fact that classically there are no long range inter-vortex forces. 

To go beyond the perturbative approximation for the Wilson and 't Hooft loops, first we ask which configurations 
are likely to give an important contribution to the path integral. The obvious guess to this are the vortex loop 
configurations. As mentioned in the introduction, with n = these loops are just Nielson-Olesen vortex loops, but 
including the CS term the form of the vortex solutions in Euclidean space is somewhat unclear. The CS term becomes 
imaginary in Euclidean space while the remaining action is real. Classical solutions will typically involve complex 
valued fields, the theory however is defined by the functional integral over real- valued fields. Our choice of path integral 
excludes any solutions to the field equations and thus any true instanton solutions. To exploit complex critical field 
configurations in the path integral requires analytic continuation of the contour in the infinite dimensional space of 
field configurations. Here we avoid this approach and restrict the path integral to real fields. If we take k to be small 
then it is reasonable to assume that the same configurations that are important for the path integral in the absence 
of k will again be important here. This amounts to evaluating the action for k^O with the vortex loops which were 
the instanton configurations solutions when k = 0. 

A similar approach was taken in [[5|, concerning 2+1 dimensional compact QED with a CS term. The authors 
assume weak coupling and treat the CS term as a perturbation. This model without the CS term was first considered 
by Polyakov, 0, who showed that the theory is linearly confining because of instantons. In this case the instantons 
which are monopoles, change the logarithmic perturbative potential between electric charges into a linear confining 
potential. Using the assumption that one may treat the CS term as a perturbation the authors of g showed that linear 



confinement no longer holds. In the instanton calculation, when the CS term is evaluated for the monopole the result 
turns out to depend on the U(l) phase of the monopole. When integrating over all possible monopole configurations 
in the path integral, one integrates over the 17(1) phases of the monopoles and these interfere destructively, cancelling 
their contribution to the path integral. In addition, since the photon is massive the perturbative contribution to 
the potential between external charges falls off exponentially with their separation. It is thus argued in that the 
confinement of electric charges does not hold when the CS term is included. In ||], the above approach of inserting 
instanton solutions from the k = case into the action for k ^ is questioned. The authors of |p| find complex 
monopole solutions with finite action and argue that that these are the relevant configurations which dominate the 
path integral, and are not cancelled by integrating over their phase variables. We feel that the treatment of the 
complex solutions in the path integral is as yet unclear since a prescription is required on how to analytically continue 
the field configuration space, and the question of gauge fixing remains problematic. Nevertheless this approach clearly 
deserves more attention. 

Complex vortex solutions are easily seen to exist in the present case, although the above controversy notwithstand- 
ing, we will not use them. We expect the situation to be less complicated than that of monopoles in the compact U(l) 
case . This is because the vortex field strength vanishes exponentially quickly outside the vortex core so we do not 
have troublesome boundary terms which could result in gauge dependence. We t ake the same ansatz for the vortex in 



Euclidean space as for that in Minkowski space considered by Paul and Khare, (2.3). With the obvious replacement 
Aq — » A' = iAo the ansatz satisfies the Euclidean equations. This is just a consequence of the fact that the vortex is 
time independent. A straight vortex line results by assuming independence on x^. The solution has finite real action 
per unit length. However the action per unit length is less than the energy of the Minkowski vortex because A' is 
imaginary. The action density is given by, 

n,.dA dA l f dA 02 e 2 22 



where A, Aq, f are the same as in (2.3). This differs from the energy of the Paul and Khare vortex in that the terms 
in the second line are negative, and that the CS term contributes to the Euclidean action, but not to the energy of 
the Minkowski vortex. The action per unit length is less than the energy of the Paul and Khare vortex by terms of 
0(k 2 ). For large k the action per unit length becomes negative. We expect here that the approximation of treating 
the CS term as a perturbation fails badly, but for small K we expect it to work well. 



III. EVALUATION OF THE WILSON AND 'T HOOFT LOOPS 



We now outline the calculation of the vortex loop contribution to the expectation values of the Wilson and 't 
Hooft loops for k =/= 0. The path integral is split into a integral over all vortex loops and the perturbative fluctuations 
from each loop. It is not possible to evaluate the exact contribution of the perturbative fluctuations, they give a 
determinant factor, which we presume is not important to our analysis. The contribution of the vortex loops to the 
path integral can be evaluated by approximating each loop as closed string whose location is given by the position 
of the zeros of the Higgs field along the vortex loop. The dominant contribution to the action of a vortex loop is its 
length times its mass per unit length. The loop may have an arbitrary shape, however if the loop twists in an unruly 
fashion or intersects itself there will be additional potential energy. This suggests in addition to the basic term of the 
mass, terms involving the curvature of the loop and a term which describes the potential energy of an intersecting loop 
could be taken into consideration. It is in principle possible to determine the detailed form of the interaction terms by 
referring back to the field theory model and examining how the equations for the Nielson-Olesen vortex change if the 
shape of the straight vortex is deformed J9| . It is also necessary to include terms describing the interactions between 
different loops. Since for a vortex loop the fields approach their vacuum expectation values exponentially quickly in 
the transverse directions, the only long range forces are those mediated by the CS term. 

It is well known that the functional integral over a gas of loops corresponds to a field theoretic path integral, p0[ , 
Jll| , H. This equivalence is demonstrated and well studied in the above papers, we will not review it here, but just 
adapt it for our present needs, where, especially, k^O. For the contribution of the vortex loops to the partition 
function e~ SB , one gets the following field theoretic path integral, 



C = N JV X V X * exp{- [[d iX diX* + mlxx* + V( X X*)]d 3 x} , (3.1) 

with J\f an normalisation factor and x is the charged vortex field. This corresponds to an effective Lagrangian for 
the vortex field x which can be treated semi-classically. The mass term for the vortex field, mj, has two contributing 
factors, one from the mass of the Nielson-Olesen vortex and an opposing factor due to the entropy factor of all possible 
vortex loops. Working on a fixed cubic lattice, for example plf , the combination m — In /i naturally emerges as being 
proportional to the mass m§, where m is the vortex mass per unit length and, a ss 5, comes from the number of 
available directions the vortex path can take. A phase transition will occur at to = In /i, where loops of all lengths 
become important. A similar transition will occur here for some value of the parameters. Once the vortex mass 
becomes small enough loops of all lengths are important and the vacuum becomes full of vortex/anti- vortex pairs. 
This is the spaghetti vacuum scenario discussed in |12| . 

The exact form of the potential V(xX*) depends on the vortex interactions. For example, repulsive delta functions 
interactions between the vortex loops which prohibit the loops from intersecting, yields a MxX*) 2 potential. To 
take into account all the correct interaction terms it is necessary to consider the vortex equations. Vector forces 
between different part of the vortex loop, and different loops, necessitates the introduction of auxiliary vector fields 
which couple in a gauge invariant manner to the vortex scalar field, x- The form of the field strength depends on 
the interaction, for example QED3 arises from Biot-Savart forces between different elements of the vortex loop. The 
vector fields, and other possible scalar fields, are then integrated out to yield the vortex potential V(xx*) which is in 
general nonlocal. 

For non-zero k most of the above discussion goes through unaltered. We saw previously that the vortex solutions 
still exist when n jt in Minkowski space, in Euclidean space any solutions are necessarily complex. We insert 
the Euclidean vortex configurations found at n = into the action when k ^ 0, which should be a reasonable 
approximation for k, small. We can assume that all vortex loops arise from vortex/anti-vortex pairs with flux 2n/e, 
since loops containing a multiple of this flux are exponentially suppressed. It is not difficult to see that evaluating the 
CS term for vortex loop configurations measures the linking numbers between different loops. For two non-intersecting 
vortex loops evaluating in/4 f cPxeijkAiFjk, gives to leading order iir 2 nn/e 2 where n is their linking number. So the 
effect of inserting Nielson-Olesen vortex loops into the complex action is to multiply their action by a phase which 
measures how the vortex loops are inter-linked. A local expression for n(C, C") is given by 



_j_ f dx f (x - y) x dy 
4tt J c ^ J c , |x-y| 



n(C,C') = -^-l dx-/ '* :/'* uy . (3.2) 



The linking number is thus manifested as a vector interaction between the vortex loops. Combined with the previous 
vector interactions from above, this will yield a vector field Gi coupling in a gauge invariant manner to the vortex field 
X- The field strength will now have a complex part, due to the fact that the linking number interaction is imaginary. 
For example, combining the linking number interaction with the Biot-Savart interaction yields a complicated nonlocal 
field strength whose lowest order terms in k are, 

1 



2i H2+j 4RF H ' VxH r-- ' (3 - 3) 

where H = V x G, and e' is a parameter which measures the strength of the Biot-Savart interaction. Again the 



correct procedure is to integrate out the vector field Gi, yielding (B.l). It is important that V(xx*) depends only on 
XX* 1 an d n °t on X or X* separately. This is due to the global symmetry x ~> & X which corresponds to conserved 
vortex number. The results below are not sensitive to the exact form of V(xx*)- Th e me chanism of how the spaghetti 



vacuum occurs is unchanged relative to the k — case. When the parameter to,q in (3.1) becomes negative the vortex 
field will acquire a nonzero vev. The vortex loops are stabilised by terms of the form X(xx*) 2 m ^(XX*) arising from 

the repulsive forces. We denote the critical value of rj 2 (as to, the vortex mass), where the phase transition occurs as 

2 

Icrit ' 

To evaluate < W(C) >—< exp(iq § c Aidxi) >, (here A, is the gauge field of the original theory), note that 
§ c Aidxi measures 2ir j e times the linking number of a vortex loop configuration, n(C), with the curve C and has the 
integral expression, 

n(C) = j> dXiPi(x) , (3.4) 



where 






jk{x - y) 1 dy k 
\x-y\ 3 



(3.5) 



The sum over C represents the sum over the vortex loops in a particular field configuration. We emphasise that this 
linking number does not come from the CS term. Thus to calculate the expectation value we do the same sum as 
before but now include the linking number factor. This essentially gives the complex scalar field theory in the external 
gauge field Pi which the Wilson loop generates, 0, O]. We get 



A/" 



J V X V X *exp{~ J[D lX D lX * + ml XX * + V( XX *)]d 3 x} , (3.6) 



where Di = d^ — 2niq/ePi. 



Note that there is no field strength for Pi, it is just an external field. (3J3) is to be evaluated semi-classically. 
When m 2 , > 0, < X >= is expected for the vacuum and the topological symmetry is unbroken. X = is a solution 
to the equations of motion even for non-zero Pi and the Wilson loop is of order 1. Vortex loops do not contribute 
appreciably to the Wilson loop. When m 2 , < however, X wu l have a non-zero vacuum expectation value, which is 
called topological symmetry breaking < X >= Xo- Inserting this as a trial solution into the path integral gives 

< W(C) >=< ex P [iq I A-dx] >« exp[-{^) 2 xl W— )] , (3.7) 

Jc e r 

where again C is rectangular with lengths r and t, t 3> r, with r representing the separation of the electric charges, 
and ro is a cut off of the order of the vortex width. For q > e non-constant screening solutions provide a better 
minimisation of the action g. The upshot is that q/e in (3.7) is replaced by Aq = q/e — n where n is the nearest 
integer to q/e. The Wilson loop has logarithmic dependence on the charge separation and is periodic in the external 
charge q with period e, which evinces the fact that when the spaghetti vacuum forms we are back in the region of 
unbroken U(l) or symmetry unbroken phase: there exist elementary charged particles of charge e which can screen 
external particles of charge e or an integer multiple thereof. This is not a conclusion that one would have expected. 
With a massive photon one would not expect any form of confinement, as is the case for compact U(l) with a CS term 
H, pi. The essential difference between the two cases is that is that when evaluating the CS term for the monopole, 
it explicitly depends on the monopole phase, and integrating over these phase variables the monopole contribution 
cancels. In the present case, the effect of evaluating the CS term for a configuration of vortex loops gives a phase 
proportional to the linking numbers of the different loops. Given a gas of such vortex loops there is no reason to 
expect these phases cancel. This is manifested by the fact that V( X X*) depends only on X *X- 

For large values of rf the vortex loops only give a small contribution to the path integral due to their large mass. As 
rj 2 is decreased vortex loops become more important until there is a phase transition at rj 2 rit where loops of all length 
become important and the U(l) symmetry is restored. For r/ 2 > i] 2 rit the Wilson loop has perimeter law behaviour. 
The perturbative contributions give a potential that falls off exponentially quickly between external charges. The 
nonperturbative vortex loop configurations may alter this potential, arising from short vortex loops which link or 
intersect the Wilson loop. But their exact contribution to the potential is undetermined from this approach since 
when the vortex loop intersects the Wilson loop there is no simple answer for the contribution to the expectation 
value. 

For 77 2 < the Nielson-Olesen vortex loop configurations do not exist. The question is whether flux tube configu- 
rations can be important when rj 2 < 0. We believe that vortex configurations continue to be relevant to the vacuum 
structure throughout the unbroken symmety phase even though they are absent classically. In [Q for the Abelian- 
Higgs model it is shown that the vortex operator possesses a non-zero expectation value throughout the Coulomb 
phase. Numerical evidence presented in Q indicates that the vortex condensate persists throughout the Coulomb 
regime. This suggests that the confining behaviour will persist throughout the unbroken symmetry phase. 

Next we turn to the evaluation of the 't Hooft loop defined along a curve C as before. Now the complex part 
of the Euclidean action plays a crucial role. To calculate the 't Hooft loop expectation value we must sum over all 
field configurations which contain an infinitely thin flux tube along C. This is equivalent to calculating the partition 
function, e~ SE , where the functional integral is over all gauge configurations A' i: such that, A\ — A. b + diVtc with 



flc as in (2.5) and Ai is non-singular. In the absence of the CS term, because the vortex loops have no long range 



interactions, the 't Hooft loop will have perimeter law dependence. Inclusion of the CS term however implies that for 
each vortex loop C" that links n(C, C") times with C we get a term, 

ex P {in(— ) 2 n(C,C")}, (3.8) 

e 

where 27r/e is the flux carried by the 't Hooft loop C. This in turn means that the 't Hooft loop calculation is identical 
to the Wilson loop calculation with the replacement q — > 2irn/e, which due to the Gauss law, is nothing else than 
the electric charge of a vortex that has flux 27r/e. Hence in the region rj 2 < rj 2 rit , vortex loops give a logarithmic 
contribution to the 't Hooft loop expectation value and one expects logarithmic confinement of external vortices. 
Such confinement of vortices has already been seen to occur in a related model fig ] . The previous remarks in the last 
paragraph about the validity of this calculation when r] 2 < also apply here. If the external vortices have electric 
charge equal to a multiple of e then the logarithmic potential disappears, however we understand this in terms of 
screening by the elementary charged scalar particles of the theory. Thus it is actually the electric charge of the vortices 
that is behind the confinement. It is again clear that we are in the region of unbroken symmetry here, since for long 
range forces to exist among the external vortices, the dynamical vortices must have disappeared from the theory 

Electrically charged particles of charge q mod e ^ are logarithmically confined in the topological symmetry break- 
ing phase (or spagetti vacuum) as discovered by Samuel Q . Our new result is that this remains true with the addition 
of the CS term, at least for small CS coefficient. Additionally, this includes the logarithmic confinement of vortices 
since they also carry electric charge due to the Gauss law. This result is surprising in view of the fact that all gauge 
fields are massive and interactions short ranged and demonstrates the importance of vortex configurations for the 
phase structure of the theory. As a note of caution, we stress that the methods employed here are far from rigorous. 
The procedure of obtaining an effective vortex Lagrangian is more problematic than when considering, for example, 
monopoles or other point-like instanton configurations, with a well defined set of zero modes. Obtaining the vortex 
action here, by summing over all string configurations has intrinsic ambiguities. Finally the question of complex field 
configurations and their relevance in the path integral needs to be examined in more detail. 
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